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The time-evolution equation of a one-dimensional quantum walker is exactly mapped to the three- 
dimensional Weyl equation for a zero-mass particle with spin 1/2, in which each wave number k 
of walker's wave function is mapped to a point q(fc) in the three-dimensional momentum space 
and q(fc) makes a planar orbit as k changes its value in [— 7r,7r). The integration over k providing 
the real-space wave function for a quantum walker corresponds to considering an orbital state of a 
Weyl particle, which is defined as a superposition (curvilinear integration) of the energy-momentum 
eigenstates of a free Weyl equation along the orbit. Konno's novel distribution function of quantum- 
walker's pseudo-velocities in the long-time limit is fully controlled by the shape of the orbit and how 
the orbit is embedded in the three-dimensional momentum space. The family of orbital states can 
be regarded as a geometrical representation of the unitary group U(2) and the present study will 
propose a new group-theoretical point of view for quantum-walk problems. 
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I. INTRODUCTION 

Quantum random walk models 0- [1 have been 
intensively studied as fundamcntal models in the bàsic 
theory of physics to discuss the relationship between de- 
terministic time-evolution with the probability interpre- 
tation in quantum mechanics and stochastic processes in 
statistical mechanics, and in the quantum information 
theory to invent new algorithms for quantum comput- 
ers (see, 0, 0, 

for recent reviews). The interest in 
novel properties of quantum walks is now spreading to 
many research fields and new applications are being re- 
ported. For example, in probability theory a new type of 
limit theorems quite different from the usual Gaussian- 
type central limit theorems was proved d @, [kJ, and 
in the solid-state physics of strongly correlated electron 
systems the Landau-Zener transition dynamics was re- 
lated to a quantum walk [l·l"j | . This paper presents a new 
aspect of quantum walks, by showing that at each time 
í € [0, oo) the quantum state of a one-dimensional quan- 
tum walker is exactly mapped to that of a Weyl parti- 
cle (zero-mass particle with spin 1/2), which is obtaincd 
by a curvilinear integration of energy-momentum eigen- 
states of a free Weyl equation along a planar orbit appro- 
priately embedded in the three-dimensional momentum 
space. In a plane polar coordinate on an orbital plane, 
(ç,7), q € [0,oo), 7 £ [— 7r,7r), the equation of orbit is 
given by 
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where a 6 C with \a\ E (0, 1) is onc of the parameters 
of the unitary matrix in U(2), which specifies the time- 
evolution of quantum walker. 

In order to show the fact that the quantum random 
walk model can be naturally considered to be a quantum- 
mechanical generalization of usual random walk mod- 
els, and also in order to clarify key points of quantum- 
walk problems, we first formulate the one-dimensional 
simple symmetric random walk as a special case of the 
following classical random-turn model. (Such a model 
is also called the correlated random walk [T^O Con- 
sider a walker at the origin of a one-dimensional lattice 
Z = {• • • , —2, —1, 0, 1, 2, • • •}, who is directed to the left 
with probability q and to the right with probability 1 —q. 
The walker has a coin, which is randomly tossed giving 
heads with probability p and tails with probability 1 — p. 
He tosses the coin and, if the outcome is head, he changes 
his direction, from the left to the right or from the right 
to the left, and if the outcome is tail, he does not changc 
his direction. In any case, then he makes a forward step. 
At the new position, he tosses the coin again, does or 
does not change his direction following the outcome of 
the coin, and then make a forward step. We assume that 
the walker repeats such random turns and steps n times 
and the probability that he arrives at the position x € Z 
and also that he is directed to the left (resp. right) after 

the n-th step is denoted by Pn\x) (resp. P„ \x)). A 
simple application of the Fourier analysis gives 



with 
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where i = — li the coin is fair (p = 1/2), the eigen- 
values of the transition matrix W(k) are A = and A = 
cosfc, and P„ (x) = P„ 2 \x) = P n (x)/2 is independent 
of q (i.e. independent of the initial direction of walker). 
This symmetric case realizes the simple symmetric ran- 
dom walks, since P n (x) — J^(dk /2n)e tkx cos n k — 

{{ n +x)/2) ■ ^et temporal and spatial units be r and 
a, respectively, and set n — t/r,x ^ x/a,k ^ ak. Then 
we consider the asymptotic of the probability density 
Pt{x) — P t / T (x/a)/a in the continuum limit a,T — ► 0. 



In the so-called diffusion scaling limit r 



0, for 



(cosafc) í/r = (1 - a 2 k 2 /2 + ■ ■ -) í/t -> e~ tk 'l 2 , wc fmd 
convergence of the probability density to 



Pt{x) 



dk 
2^ 



exp 



t 



-k + ikx 



-x 2 /2t 



which is called the heat-kernel, since it sol ves the heat 
equation [d/dt — (l/2)d 2 /dx 2 ]p t (x) = with the initial 
condition lim t _,o Pt{x) = 5(x). The limit of walker 's po- 
sition is in the Gaussian distribution with mean zero and 
variance t, and this convergence property of random-walk 
distribution is a typical example of general central limit 
theorems. 

Now we introduce the quantum random walk model as 
a quantum version of the above random-turn model. In 
the LHS of (J2J) , the two component vector is replaced by 
a two-component wave function 



^n(x) 



* { n\x) 



(4) 



and in the RHS of this equation, the initial distribution 



of walker's directions 



by an initial qubit 



1-qJ y p 
a,(3 G C, \a\ 2 + \f3\ 2 = 1, the transition-probability ma- 
trix W(k) by a transition matrix of probability ampli- 
tudes 

ik 



U(k) = 











A. 



(5) 



Here A is a 2 x 2 unitary matrix, which determines be- 
havior of a "quantum coin" . 

In the context of quantum mechanics, the two- 
component wave function usually describes quan- 
tum state of a particle with spin 1/2, and the Hilbert 
space of spin operators is spanned by the Pauli matri- 
ces <Jj,j — 1,2,3, and the unit matrix I2 of size 2. In 
the classical random-turn model, the direction of step 
is determined by the walker's direction that was ran- 
domly changed just before doing the step. Correspond- 
ing to this setting, operators which determine motions of 
quantum walkers should be coupled with spin operators 
acting "spin states" of the walkers. In the usual three- 
dimensional quantum mechanics, the momentum vector 
operator p = (pi,P2,P3) = —iK^ is the operator for the 
particle motions and the easiest coupling with spin oper- 
ators is given in the form 



where er = (ci, (12,(73). This is known as the Hamiltonian 
for free motions of a zero-mass particle with spin 1/2 like 
a neutrino. The corresponding Schròdinger equation 



d , 



H(p) 



er p, 



(6) 



Íh—$ t (p) = (cr ·p)$ t (p) 



is called the Weyl equation in the three-dimensional mo- 
mentum space [13f , which is reduced from the Dirac equa- 
tion for four-component wave functions by setting the 
mass term to zero (see, for example, Section 10.12 of 01 
and Section 2-4-3 of 

In the present paper, we will show that the time- 
evolution equation of a quantum walker in the fc-space 
(one-dimensional Fourier space) is exactly mapped to the 
Weyl equation in the p-space (three-dimensional momen- 
tum space). Since the quantum random walk model is 
defined on a lattice Z, the wave number k should be in a 
finite region [— tt, tt) with periodicity k±2ir = k (the Bril- 
louin zone). In other words, the fc-space can be identificd 
with a unit circle, in which k describes an angular coordi- 
nate of a position on the circle. We will show that, by our 
nonlinear mapping fc i— > p(fc) = (pi(k),P2(k),pa(k)), this 
unit circle is transformed into a planar orbit in the three- 
dimensional p-space. The obtained orbit is no longer a 
unit circle, and the radial coordinate q depends on the 
direction on the orbital planc. If wc dcfine the angular co- 
ordinate 7 G [— 7r, tt) by an angle from the direction to the 
"perihelion" of orbit from the origin, the equation of orbit 
is given by The deformation of the unit circle, with 
which it is embedded in the p-space, is fully described by 
the Jacobian J = \dk/d~y\, which appears when we con- 
sider the transformation of the integral f* 7r (dk/27r)f(k) 
to the curvilinear integral along the orbit in the p-space. 
One of the main results reported in the present paper 
is that Konno's function 0, , which describes the dis- 
tribution of the pseudo-velocities of quantum walker in 
the long-time limit n —* 00, is directly related with this 
Jacobian. Moreover, it will be shown that, the normal 
direction of the orbital plane depends on the choice of 
unitary matrix A, and from this dependence the initial 
qubit dependence is completely determined. 

This paper is organized as follows. In Sec. II, the precise 
description of the quantum random walk models and the 
bàsic formulae for physical quantities studied in this pa- 
per are given. Konno's weak limit theorem is then briefly 
reviewed. We rewrite U(k) by using the exponential op- 
erators with the Pauli matrices, and then the exact map 
to the Weyl equation is given in Sec. III. There the one-to- 
one correspondence between the quantum-walker states 
in the fc-space, each of which is specificd by the unitary 
matrix A, and the orbital states of a Weyl particle in the 
three-dimensional p-space is clarified. Then in Sec. IV we 
follow the argument by Grimmett et al [Toj and give a 
new proof of Konno's weak limit theorem. Concluding 
remarks are given in Sec.V, and Appendices are used for 
some details of calculations. 
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II. THE MODEL AND KONNO'S 
DISTRIBUTION 

A. Quantum random walk models associated with 
unitary matrices 

Consider a two-component wave function Q). Follow- 
ing the definition of the model given by Konno 0, , we 
have 

* ( r )- (^n\x + l)+b¥^(x+l)\ ( ) 



(i b 

c d 

matrices. If we set 



where A = ( " ) £ U(2) e a set of 2 x 2 unitary 



and assume the relations 

*«(fe) = £*<?>(x)e- ite , j = l,2, 

then O is rewritten in the wave-number space (fe-spacc) 
of the Fourier transformation, fe £ [— n, tt) as 

*„+i(*) = U(k)Ú n (k), ti = 0,1,2,- ••, 

where U(k) is given by JSJ. 

The state at time-step n is then obtained from the 

initial state * (fe) = í " J , a, /3 € C, |a| 2 + |/3| 2 = 1, by 



*„(fc) = C/(fc)"* (fc), 



(8) 



whose Fourier transformation gives the wave function at 
time-step n as 

í/k 

g^^fc) 
^ e ^E/(fc)«v£r (fc). 

The distribution function in the real space at time-step 
n is given by 

Pn(x) = |*„(z)| 2 = *+(a:)* n (a;) 



c i(k-k')x 

2tt 2tt 



— 7T «/ — 7T 



^(feO^*O n )(W*o(fc)). (9) 



It is known that any matrix in U(2) can be parameter- 
ized as a matrix in SU(2) times a phase factor in the form 
e iv , tp £ [-vr/2, 7t/2), where SU(2) denotes the set of 2 x 2 



unitary matrices with determinant 1. For example, the 
Hadamard matrix 



F = 7i0 -i) eu(2) 



(10) 



is written as H = e^^A with A = -4 I ".] 6 



75 \ i -i 

SU(2). Since the phase factor e lip of U(k) is irrelevant 
in calculating the distribution function JÏÏJ), without loss 
of generality, we can assume that the matrix A in is 
chosen from SU(2). We can see that SU(2) is generally 
written as 



SU(2) 



1=1 l );a,òeC,|a| 2 + |6| 2 = 1 



= iA = 



a 



- u 2 e l4 

,-i0 



u £ 



[O,l],0,^e[-7r,7r)}. (11) 



As shown by the second equality above, it is a three- 
dimensional space parameterized by real variables u, 9 
and <j) (Cayley-Klein parameters). 



B. Formulae of expectations 

Let X n denote the position of the one-dimensional 
quantum walk at time-step n — 0, 1, 2, • • ■ and consider a 
function / of x £ Z. The expectation of f(X n ) is defined 
by 



/(*«)) = X;/(aOP„(aO 

x£Z 

= £/(*) f f^-^(fe') f ^e«-* n (*). 
If f(x) = x r , r = 0, 1, 2, • • •, it is written as 



x£Z 1 

{(-4) «"■>*•<*>■ 



We note that $ n (fe) should be a periòdic function of fe £ 
[— 7T, 7r), and then by partial integrations, we will have 



and thus 
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where we have used the summation formula X^ez e = 
2ir5(k). Then, if f(x) is analytic around x = 0, triat 
is, if it has a converging Taylor expansion in the form 



with 



m = e ■ 



c J , we will have the formula 



r/A; 



f(X n ) = — i— *„(*) 



2tt 



cZ/r 



(12) 



C. Konno's results 

A fundamental requirement of the quantum mechanics 
is that any time-evolution of an isolated system is given 
by a unitary transformation of wave function, which is 
necessary to allow us to adopt the usual formula for prob- 
ability distribution of outcomes in observations given by 
using squares of the wave function. In the present case, 
U(k) is unitary and the probability that the quantum 
walker is observed at the position x £ Z after n-th step 
is given by (JÜJ. The expectation of a function / of posi- 
tion at time-step n is then calculated as <|12|) . 

It should be noted here that the unitarity of U (k) en- 
sures |A| = 1 for any eigenvalue A of it. This fact implics 
that in principle we are not able to find any convergence 
property of wave function ^f n (x) nor of the probability 
P n (x) in the long-time limit n — > oo in quantum- walk 
problems. It presents a striking contrast to the classical 
systems of random walks, which will generally converge 
to diffusion particle systems in the long-time and large- 
scale limit, as we have demonstrated in Sec.I in the simple 
and symmetric case. 

Recently, however, Konno proved that the distribution 
of the pseudo-velocities of the quantum walker X n /n does 
converge in such a weak sense that any moment of the 
pseudo-velocity converges in the long-time limit n — > oo 
to a moment of a random variable, whose distributed is 
given by a novel probability density function |g, |9( ■ That 
is, for any r = 0, 1, 2, • • •, 



((Xn/r 



dy y r v(y) 



where 



\a\) = - 

""(i - y 2 )v \ a \ 2 - y 2 

I(y;a,b;a,(3) 



(14) 



V- (15) 



Here denotes the indicator function of a condition 
lo; = 1 if ui is satisfied and = 0, otherwise. The 
special case, in which the matrix A is (proportional to) 
the Hadamard matrix Q10[l . has been studied well with 
the name of Hadamard walk, since it corresponds to the 
p = 1/2 case of the classical random-turn model @ with 
(0. But Konno's weak limit theorem claims that even in 
this case, 1^1; the limit distribution depends on the 
initial qubit and it is not symmetric in general. The limit 
distribution is very sensitive to the initial qubit, and if 
and only if the initial qubit is chosen as |a| = and 
af3*ab* + a* (3a* b = for given a and b in A, probability 
density becomes a symmetric function n(y; |a|)l{| a |<| a |} 
(see[|i). 



III. MAP TO WEYL EQUATION 
A. Exponential operator with Pauli matrices 

Consider the Pauli matrices, o\ = ( ^ g j , <72 



For an arbitrary three- 



-i\ (l 

1 o y aa - [o -i, 

dimensional vector q = (?i, 92, 93), we can see that 
(er • q) 2 = q 2 h, where q = |q|. Then 



00 

= £ (- í<T . q )« 

71 = 

= \ I2 — icr ■ q tang cosg, (16) 



v(y) = /j,(y; \a\)T(y; a, ò; a, /5) 1 -[ | < | o.| } (13) where q is a unit vector defined by q = q/g. 



By using the well-known àlgebra of the Pauli matrices, we find that JSJ) with A £ SU(2) can be written as 



U(k) = 



11 e 



ucos(Zc + 9) 



.. \/l-u 2 sin(Zc + 4>) \J\-u 2 cos(Zc + (f>) 
h + i\ ; ; . „ ; + ; ), , + tan(fc + 0) a 3 



u cos(fc + 0) 



u cos(k + 6) 



(17) 



In order to identify U(k) with (|16|l by choosing the vec- ucos(fc + 6) = cosq, (vï — u 2 /u)(s'm(k + </>)/ cos(fc + 
tor q appropriately, we consider a system of equations; 
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9)) = -^tang, (-y/ï - u 2 /u)(cos(k + 0)/cos(fc + 6»)) 
— <?2 tanç, tan(fe + 0) = —(73 tang. It is solved as 



q(k) 



arccos 



= arctan 



wcos(fc + 9) 
1 



cos(fc + 9) V u 2 



— - cos 2 (fc + 9) 



(18) 
,(19) 



where we take the principal value of arccos x and choosc 
appropriate branches of arctan x so that q{k) is a periòdic 
and continuous function of k G [— 7T, 7r), and 



&(*) 



Vi 


- u 2 /it 


VV« 2 - 


cos 2 (fc + 9) 


Vi- 


- v? ju 


Vl/u 2 - 


cos 2 (fe + 9) 




1 


y/l/u 2 - 


cos 2 (fe + 9) 



sin(fc + (/>), 
cos(fc + (/)), 
sir±(fc + 0). (20) 



We then define a threc-dimcnsional vector 

q(fe) = (a{k)qí{k),q{k)q2{k),q{k)q^k) 
By the formula (|16|l . [/(fe) is now written as 



(21) 

(22) 
to define 



U{k) = e-^'iW. 

We can interpolate time-steps n = 0,1,2,' 
a continuous time t G [0, 00), and we have 

* t (fc) = e- iíO "· q W* (fc), fe G [-7T, tt), (23) 
which solves the Weyl equation with % = 1, 

^*í(fc) = W(q(fc))*t(fc), (24) 

where the Hamiltonian is given by The equations 
(|18[l - <|21[l define a map k G [— 7r, 7r) i—* q G p-space, and 
by this map the time-evolution equation © of the one- 
dimensional quantum walk is exactly transformed to the 
Weyl equation . 

It is easy to solve the eigenvalue problem of the Hamil- 
tonian (JüJ) of a free Weyl particle with momentum p 
[Ï4L HE) . The eigenvalues are A = ±p, where p = |p|. 
The corresponding energy-momentum eigenfunctions are 
independent of the absolute value of p and given as func- 
tions of p = p/p by 



VW(p)/ 



V'-(p) 



>«(p)' 
^ 2) (p) 



■ sin(6» p /2)e-^" 
cos(0 p /2) 



(25) 



in the polar coordinate of momentum, p\ = 
p sin 9 p cos ip p , P2 = p sin 9 p sin ip p , p^ = p cos # p , which 



are orthonormal in the sense that 

I^+(P)| 2 = ^(P)^+(P) = 1, 

|VMp)| 2 = vA(p)?Mp) = i, 

^(P)^-(P) = V'1(P)V'+(P) = 0. (26) 

Now the problem is reduced to the study of the property 
of nonlinear map k n- q G p-space. 

B. Planar orbital states in the momentum space 

Since (|18fl - i|20|) show that all the components qj (fe) , j = 
1,2,3 are periòdic functions of k G [— 7r, 7r), q(fc) gives 
a one-parameter orbit in the three-dimcnsional p-space 
parameterized by k G [— tt,tt). Define a unit vector 

è 3 = (-u cos(0 -9),u sin(0 - 0), y/ï - w 2 ) . (27) 

Then we can see q(fe) • §3 = 0. That is, 

q(fe) _L §3 for all fcG [— 7r,7r), 

which implies that the orbit is on a plane including the 
origin, whose normal vector is è^. We denote this orbital 
plane by H(u, 9, (f). 

Since we have assumed the principal value for arccos x, 
< arccos u < tt/2 for any given < u < 1, where 
arecosw = 7r/2 if u = 0, and arecosw = if u = 1. As 
explained in Appendix A, we can see that q(k) — |q(fc)| 
takes its minimum value g m i n = arccos u when k = —9 
(mod 27r) and its maximum value (7 max = i" — <7min when 
fe = -tt - 6 (mod 2tt). When fe = -tt/2 - 6» and k = 
tt/2 - 9 (mod 2tt), g(fc) = tt/2. Let 



ei 

é 2 



0),-cos(<£-0),g), 



v/r 



),-vT 



(28) 



so that a set (èi, 02,63) makes an orthonormal basis in 
the p-space and the orbital plane H(u, 0, <p) is equal to 
the (èi, è2)-plane. The unit vector èi is directed to the 
"perihelion" q(feo) with fco = — & (mod 2tt) and we define 
the angle 7 in this plane by 

cos 7 = q(fc) • èj., 

where q(fc) = q(k)/q(k). Then we have 



cos 7 



S1117 



(\Jl-u 2 /u) cos(fc + 0) 
y/l/u 2 - cos 2 (fc + 0) 

(1/u) sin(fc + 0) 
x/l/u 2 -cos 2 (fc + 0)' 



(29) 
(30) 



Combining i|19|) and (|29|l . we have the equation for the 
orbit 



tanç 



y/1 - u 2 1 
u cos 7 



(31) 
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in the plane polar coordinates (q, j),0 < q < 00,7 S 
[— 7r, 7r), on the orbital plane H(u,9, </>). We can rewrite 
(I18|l - H21ll as the functions of q and 7 as 



91 = gvi- 

g 3 = qu sin7 



COS — 



) SIE7 — q sm(ç> 
9) shi7 — gcos( 



) cos 7, 
0) cos 7, 
(32) 



which are the representations for the three components 
of q by the plane polar coordinates (q, 7) on H(u, 9, 4>). 

It is interesting to see dependence of the orbit 
{q(fc);fc G [— 7r,7r)} on the parameters u,6,<p of uni- 
tary matrices. The normal direction of orbital plane 
63 given by (|27() is in the z-direction, when u = 0, 
and it is tilted, as u increases from to 1, to the 
(— cos(0 — 9), sm(4> — 9), 0)-direction lying in the (x, y)- 
plane. We assume that arctanx takes a principal value 
(0 < arctanx < 7r/2) for x > 0, and a value in (7r/2,7r) 
for x < 0. Then H31JI gives 



q = arctan 



vT 



1 



cos 7 



When u = 0, q = 7r/2, which implies that the orbit is a 
circle with the radius 7r/2 and the center at the origin. 
Since COS7 > (resp. COS7 < 0) for 7 £ [— 7r/2,7r/2] 
(resp. 7 G [— 7T, — 7r/2) U (tt/2,tt)), in the limit u — * 1, 
we will see that q — * when — 7r/2 < 7 < 7r/2 and 
<7 — > 7r when — 7r < 7 < — 7r/2 and 7r/2 < 7 < 7r. When 
u G (0, 1), the orbit is a deformed circle. Therefore, if we 
observe the orbit riding on the tilting orbital plane, its 
shape is changing in u as shown in Fig. ^ 




=0.001 

=Hadamard 
=0.999 



FIG. 1: (Color online) The dependence on u of shape of the 
orbit. Three cases u ~ 0, u — \/\/2 (the Hadamard matrix), 
and u ~ 1 are shown. 



C. Curvilinear integration along the orbit 

Now we consider the integration with respect to fe on 
-7r, 7r), which is necessary to calculate the expectations 



(|12|l . The above result implies that this integration cor- 
responds to the integration over 7 G [— tt, tt), which per- 
forms the curvilinear integration along the orbit. We 
write the Jacobian associated with the map fe — + 7 as 
,/ = \dk/dry\. By simple calculation, we have (see Ap- 
pendix A), 



J = 



VT 



1 — u 2 sin 7 



and the relation 



2t7 j y ' J_„ 2n 1 — u 2 sin 2 7 



(33) 



/(%)) (34) 



is established, where fe (7) is determined by l|A2|) and 
(EU). 

Then we change the variable as 



to have the equality 



y = u sui 7, 



(35) 



* d ±f(k)-2 f * 



1 



/(*(!/))• 



-y 2 1-2/ 
That is, we arrive at the formula 

-/(fe)= / dy /*(»;«)/(%)) 

with Konno's function l|14|l . 

IV. LIMIT THEOREM OF QUANTUM WALKS 

A. Quantum walk starting from the origin 

We assume that the quantum walk starts from x = 0, 
that is 

T , r„\ _ xr„\ ( a \ ^ ,t, m _ Z' " 





*„(a0 = S(x) I " 1 ^ *o(fc) = 



where a,/3 G C with |a| 2 + |/3| 2 = 1. We can express 
this initial wave function í'o(fe) a s a linear combination 
of the two eigenfunctions (|25|l of the Hamiltonian H(p) 
with p = q(fc) in the way 

* (fc) = C + (q(fc))V + (q(fc)) + C_(q(fc))^-(q(fe)), (36) 

where C+(p) and C_(p) are given as functions of p = 
P/|p| by (jBlf) in Appendix B. Then, for n = 0, 1, 2, • • •, 
(|23j) gives 

#„(fe) = e -iW(q(*))f> 

x | C+ (q(fc))^ + (q(fc)) + C_ (q(fe))V- (q(fe))| 

= e-^ fc )"C + (q(fe))^ + (q(fe)) 

+ e ^( fe )"C_(q(fc))^(q(fc)), (37) 
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since ip±(p) are eigenfunctions of H(p) with the eigen- B. Weak limit theorem 

vàlues A = ±p. It should be noted that the timc 
dependence and the initial-qubit dependence are sepa- 

rately controlled by the absolute value of the vector q(fc), We recali the formula l|12l) . It is enough to considcr 

q(k) = |q(fc)|, and by its direction, q(fc), respectively. the case f(x) — x r 7 r = 0, 1, 2, • • • . 



Following the argument by Grimmett et al. ^(j> we fi n d from l|37|l . 



~dk ) " w \ dk ) 

■^jP) e ï «W"C-(q(fc))^-(q(fc))n r + 0(n ï - 1 ). (38) 

Insert (I37|l and (|38|) in the formula l|12fl . use the orthonormality (|26p. and take the limit n — > oo. Then we have [ïo| . 

hm ((X„/nD = f ^{|a + (q(fc))| 2 + (-ir|C_(q(fc))|4 f . . = 2 ±ft ^ , (39) 
»^°° 2tt ^ J y ^1/u 2 - cos 2 (fc + 6») y 

where we used (|A1|) . 

Now we perform the map fc — > 7. By l|30l) . 1341) and (|B2|) . we have 



liii'·.l.V,,/».^: = r^ y\^ 2 (usin 7 ) 2 ™ 

27r 1 — ü sm 7 



hm ((X„/n) 2m+1 ) = - J (|a| 2 - |/3| 2 ) + ^LjZ( a/3 * e -^-<>) + a */3e^- e )) 



T c? 7 Vl - u 2 . , 2m+2 
/_„■ Z7r 1 — w z sm 7 

for m = 0, 1, 2, • ■ •. Then we change the variable as 13 5|) . The rcsults are expressed using i|14|) as 



hm ((X„/n)^) = / dy^y;u)y 



hm ((X n /n) 2m+1 ) = - { (\a\ 2 - \(3\ 2 ) + * (ap* e-*^ + a* pe*^) / dy u)y 2m+2 . (40) 



Then we have arrived at the following limit theorem, 
which was first obtained by Konno 0,Í|. 

Theorem Lct X n be the one-dimensional quantum 
walk associated with the matrix 

with a, b £ C, \a\ 2 + \b\ 2 = 1, which starts from the state 

Mx) = S(x)(j j Y \a\ 2 + \(3\ 2 = 1, a,P£C. 
Then for any analytic function f(x) on Z, 

/oo 
dyf(y)v{y) in n -> 00, 
-OO 

where z^(y) is given by (|13|I - H15|I . 



V. CONCLUDING REMARKS 

We remark that the quantity h — cr ■ p is called the /ie- 
licity. The energy-momentum eigenfunctions ^+(p) and 
tp- (p) are also the eigenfunctions of the helicity with the 
eigenvalues +1 and —1, respectively. The difference be- 
tween the behaviors of even moments and odd moments 
in 14U|I comes from the fact that these two states with 
opposite helicity contribute with the same sign to even 
moments and with alternating signs to odd moments, re- 
spectively, as we can see in (|39|) . In the quantum spin 
systems, the term of the form cr ■ H and cr ■ d appear 
in Hamiltonian, representing the Zeeman interaction be- 
tween a spin and an external magnètic field H, and the 
spin-orbit interaction with an anisotropy vector d, re- 
spectively. We expect that these analogue with other 
physical systems will be useful, when we consider many- 
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body systems of quantum walkers in the future. 

The present one-dimensional quantum walks make a 
family of models, each of which is specified by a 2 x 2 uni- 
tary matrix. The corresponding family of orbital states 
in the three-dimensional momentum space studied in this 
paper can be regarded as a geometrical representation of 
the unitary group U(2). This point of view will provide a 
useful sight in quantum-walk problems not only in solv- 
ing problems but also in setting problems themselves and 
in classifying results. For example, what kind of quantum 
random walk models can be associated with the unitary 



J 



group SU(3), in which the so-called Gell-Mann's eight 
matrices (see, for example, |ïq]) play the similar role to 
the Pauli matrices in SU(2) ? In the recent papers, a lot 
of interesting phcnomena have been reported in a variety 
of quantum random walk models, which are associated 
with 3 x 3 or larger matrices |171 E 03. Symmetry 
and invariance, which are sometimes hidden in a vari- 
ety of phenomena observed in quantum walks, should be 
clarified using the group-theoretical investigation in the 
future. 



APPENDIX A: ON THE ORBITS 



Here we assume that k + 9 e [— tt, tt) for simplicity of description. By i|18fl . we can see that 

TT TT 

(i) q(k) = tt — arccosu \ — \ arccosu as k + 9 = — tt / — — /* 0, 

TT TT 

(ii) q(k) — arccosw / — f tt — arccosu as k + 6 = y — / tt, 

where x = a /" 6 (resp. x = a \ 6) means the variable x monotonically increases (resp. decreases) from a to b. That 
is, dq(k)/dk < 0, when — tt < k + 8 < 0, and dq(k)/dk > 0, when < k + 6 < —tt. By this observation and the general 
formula (d/dx) arccosx = =pl/Vl — x 2 , we can determine the derivative as 

dq[k) sin(/c + 9) 



dk 

The equations (HHj) - (f201> g ive 

q(-7T - 0) = <7max(sil 



\Jl/u 2 — cos 2 (k + 1 



(Al) 



q(-7r/2-0) = | (yï - u 2 cos(0 - 6»), - y/í - v? 



sm 



),u 



qi-0) 



where q m j n = arccos u and g max = tt — q n 
The relations (23 and (j20l imply 



1, cos( 



,0 , 



q(TT/2-0) = |(-VT 



),V1 



-u , 



By differentiate the both sides of (|A3() . we have 



k + e 


= —TT 


44- 


7 = 




k + e = 


-tt/2 




7 = 


tt/2, 




9 = 




7 = 


o, 


k + e 


= tt/2 


4» 


7 = 


-tt/2, 



and they can be inverted as 



sin(fc + 9) = - 
cos(k + 9) = 



(Vl — u 2 /u) sin7 
— u 2 )/u 2 + cos 2 7' 
(I/u) cos 7 



— u 2 )/u 2 + cos 2 7 



(A2) 
(A3) 



sin(A; + 6>)dfc = — 
(27 



Here we see 



d_ 



(1/u) cos 7 



(1/u) COS7 



— u 2 )/u 2 + cos 2 7 



v(T — u 2 )/u 2 + cos 2 7 



CÍ7. 



sin(fc + é*) 



1 — u z sm 7 



where we have used (|A2ll . Then the Jacobian J 
\dk/d-f\ is determined as (13311 . 



APPENDIX B: INITIAL STATE DEPENDENCE This can be represented as a linear combination of the 

two eigenfunctions of 7í(p) = cr ■ p given by (|25ü as 

Here we consider a two-component unit vector 

^o = C + (p)^+(p) + C_(p)^_(p), 



a p ), a,f3eC, |a| 2 + |/3| 2 = I. 



J 



where 



C+(p) = ^(p)^ = a cos^ + /?sin^e-^, 

C- (p) - VA (p)^o = -« sin ^e** + /3 cos ^ (Bl) 



by the orthonormality H2bjl . By straightforward calculation we find 



C±(p)| 2 = \ ± i J (|a| 2 - |/3| 2 )p 3 + a/3* (pi + *Pa) + a*/3(fr - zp 2 ) 



2 2 

Now we set p =>■ 4(7) , which is given by ijïï^ . Since 



Pi±ip 2 <?i ± iq 2 = ^\/ 1 — u 2 sin 7 =F i cos7je =p ^ e \ 

the above results become 



|C ± (q( 7 ))| 2 = \±\U\a\ 2 -\P\ 2 ) + ^^{aP*e-^ + a*pé 
2 2 1 u 



>u SII17 



T^H'e-*- 6 ' - a*/3e l(</, - e) ) cos 7 . (B2) 
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